On the two-dimensional state in driven magnetohydrodynamic turbulence 
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The dynamics of the two-dimensional (2D) state in driven tridimensional (3D) incompressible 
magnetohydrodynamic turbulence is investigated through high-resolution direct numerical simula- 
tions and in the presence of an external magnetic field at various intensities. For such a flow the 2D 
state (or slow mode) and the 3D modes correspond respectively to spectral fluctuations in the plan 
fen = and in the area fcii > 0. It is shown that if initially the 2D state is set to zero it becomes 
non negligible in few turnover times particularly when the external magnetic field is strong. The 
maintenance of a large scale driving leads to a break for the energy spectra of 3D modes; when the 
driving is stopped the previous break is removed and a decay phase emerges with alfvenic fluctua- 
tions. For a strong external magnetic field the energy at large perpendicular scales lies mainly in 
the 2D state and in all situations a pinning effect is observed at small scales. 

PACS numbers: 47.27.Jv, 47.65.-d, 52.30.Cv, 95.30.Qd 
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I. INTRODUCTION 

A variety of astrophysical plasmas is well described by 
magnetohydrodynamic (MHD) turbulence in the com- 
pressible or even simply in the incompressible case [U-Q . 
The solar wind is often cited as an example since many in 
situ data are available which have shown a medium char- 
acterized by turbulent fluctuations over a large band of 
frequency / from a fraction of mHz to kHz. Whereas the 
turbulence regime at / > 1Hz requires the introduction 
of non MHD processes like dispersion [H, H| the low fre- 
quency part at / < 1Hz is the source of a lot of activities 
on MHD BM. 



The first description proposed by Iroshnikov and 
Kraichnan (IK) [33, HH is a heuristic model for incom- 
pressible MHD turbulence a la Kolmogorov where the 
large-scale magnetic field is supposed to act on small- 
scales as a uniform magnetic field, leading to counter- 
propagating Alfven wave packets whose interactions with 
turbulent motions produce a slowdown of the nonlin- 
ear energy cascade. The typical transfer time through 
the scales is then estimated as t^/ta (instead of r„; 
for Navier-Stokes turbulence), where r„; ~ Ijug is the 
eddy turnover time at characteristic length scale I and 
ug is the associated velocity. The Alfven time is the 
time of collision between counterpropagating wave pack- 
ets and is estimated as ta ~ i/Bo where Bo represents 
the large-scale magnetic field normalized to a velocity 
(Bo — >• Hoy/nopo). (This renormalization will be used 
in the rest of the paper.) Hence, the energy spectrum in 
£-3/2 un K]j C fc- 5 /3 Kolmogorov one for neutral flows. 
The weakness of the IKs phenomenology is the apparent 
contradiction between the presence of Alfven waves and 
the absence of an external uniform magnetic field: the 
external field is supposed to be played by the large-scale 
magnetic field but its main effect - i.e. anisotropy - is 
not included in the description. 



Two fundamental evolutions in MHD turbulence have 
been made during the last two decades. There are both 
concerned with anisotropy. The first one is the conjecture 
that the refined times r n i ~ £±/u£ ± and ta ~ ^\\/Bq 
(where _L and || are respectively the perpendicular and 
parallel directions to the mean magnetic field Bo) are 



balanced at all scales 



It leads to the heuristic k 



-5/3 



; 2/3 



energy spectrum as well as the relationships fcii 
Whereas the first relation is a trivial consequence of the 
conjecture (and often wrongly interpreted as the main 
result of the conjecture) the second prediction reveals a 
non trivial character of MHD turbulence. The second 
fundamental evolution is the possibility to handle the 
effects of a strong Bq on the MHD dynamics through 
a rigorous mathematical treatment of weak turbulence 
which leads asymptotically to a set of integro-differcntial 
equations. The exact solution for weak turbulence at 
zero cross-helicity is a kj 2 energy spectrum [lij]. Note 
that the form of the energy spectrum in the regime of 
strong turbulence is still the subject of discussions 

fiod although the relationships k\\ ~ k 2 / 3 seems to be 
often verified, whereas the weak turbulence prediction 
has been obtained recently by two independent set of 
direct numerical simulations @, HU • 

The transition from strong to weak turbulence has 
been the subject of a previous paper where the regime 
of decaying turbulence has been investigated [24[ ■ In the 
present paper the focus is made on the regime of driven 
incompressible MHD turbulence under the influence of 
a uniform magnetic field at various intensity and in the 
balance case. A set of high-resolution direct numerical 
simulations in the tridimensional case are reported. Par- 
ticular attention is turned to the dynamics of the 2D state 
(also called slow mode) which corresponds by definition 
to the fluctuations in the plan k\\ = 0. In particular, 
we demonstrate the fundamental role of the 2D state at 
large B . 
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The organization of the paper is as follows. In Section 
HI! the numerical setup is given as well as the definition 
of the 3D modes versus the 2D state. In Section Hill the 
temporal characteristics of the different flows are investi- 
gated whereas the study of the energy spectra is exposed 
in Section IIVI Finally, a summary and a conclusion are 
given in the last Section. 



II. NUMERICAL SETUP 

A. MHD equations 

The incompressible MHD equations in the presence of 
a uniform magnetic field Bq read 



d t v - B d\\b 
d t b - B 9||V 



v • Vv = -VP, - 
v Vb = b • V v 

V- v = 0, 

V ■ b = 0, 



b- Vb 

- ??Ab , 



/Av . 



(1) 
(2) 
(3) 
(4) 



where v is the velocity, b the magnetic field (in velocity 
unit), P* the total (magnetic plus kinetic) pressure, v the 
viscosity and r\ the magnetic diffusivity. The introduction 
of the Elsasser fields z ± = u ± b for the fluctuations is 
also useful; it gives (assuming a unit magnetic Prandtl 
number, i.e. v = rj) 

T P 3|iz ± + z T • Vz± = -VP* 



0, 



i/V 2 z ± 



(5) 

V ■ z ± = . (6) 

The third term in the left hand side of Eq. ([5]) represents 
the nonlinear interactions between the fields, while 
the second term represents the linear Alfvenic wave prop- 
agation along the Bo field which defines the z-direction. 
Note that both descriptions Q-flD and flU)-® will be 
used through the paper. 



B. Initial conditions 

We numerically integrate the 3D incompressible MHD 
equations (U}-® m a 27r-periodic box, using a massively 
parallel pseudo-spectral code including de-aliasing and 
with a spatial resolution of 512 x 512 x 512 grid-points. 
The time marching uses an Adams-Bashforth / Cranck- 
Nicholson scheme which is a second-order finite-difference 
scheme in time [35| . A unit magnetic Prandlt number is 
chosen i.e. v = ?y for all runs. We present below the 
different set of simulations. 



1. Runs la to Ilia 

The first set of numerical simulations is character- 
ized by an external force which fixes the bidimensional 
Elsasser spectra as 



where only the two first perpendicular and parallel 
wavenumbers are energized with P(l) = F(2) = 1/2. 
The 2D state, E ± (k±,0), is never forced which means 
that initially it has no energy and may evolve freely at 
time t > 0. 
The kinetic 



E v 



1 



and magnetic 



1 



<u 2 (x) >, 



<b 2 (x) >. 



E b 



(8) 

(9) 
1/2. 



energies are initially equal such that E v 
Note that < • > means a space averaging. 

The reduced cross-helicity between the velocity and 
magnetic field fluctuations which is measured by 



_ 2 < u(x) • b(x) > 
P= < u 2 (x) +b 2 (x) > ' 



(10) 



is initially set to zero. The initial (large-scale) kinetic and 
magnetic Reynolds numbers are about 1775 for the flows 
with v = 2 x 10 -3 (see Tabic HI, with u rms = b rms = 1. 
The perpendicular integral length is 



L i = 2tt 



J(E v (k ± )/k ± )dk ± 
fE v (k ± )dk ± 

whereas the parallel integral length is 



L 



2tt 



/(£«(fc||)/fc||)rffc|| 



3.5, 



5. 



(11) 



(12) 



A parametric study is performed according to the in- 
tensity of P - Three different values are used, namely 
Bq = 1,5 and 15. All these simulations are driven up to 
time t* for which a state of fully developed turbulence 
is reached. Then, the driving is stopped and the flows 
evolve freely. These simulations correspond respectively 
to runs la to Ilia (see Table IJ). 



2. Run Illb 

The second type of simulation is characterized by an 
external force on the bidimensional Elsasser spectra such 
that 



E ± (k ± ,k ll ) = F(k {l )k 2 A 



(13) 



P ± (fc_ L ,fc||) = F(k\\)k 



3 

± ■ 



(7) 



where F is the same as in 0. The kinetic energy is 
initially larger than the magnetic energy with E v ~ 2 x 
E b = 0.67 and only the case Bq = 15 is considered. 

All information presented in this Section is summa- 
rized in Table HI Additionally, we give the computa- 
tional parameters when the stationary phase is reached 
for which a balance is obtained between forcing and dis- 
sipation. 
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TABLE I. Computational parameters are given for runs la to Ilia for which initially E v = E b and for run Illb for which 
initially E v > E b (see text). Spatial resolution, viscosity v [y = rj) and magnetic field intensity Bo are given, followed by 
the initial values of the integral length scales (perpendicular L± — 2n f (E v (k±) /k±)dk±/ f E v (k±_)dk± and parallel Lu = 
2ty J (E v (ku)/ku)dku/ f E° (kn)dku), the r.m.s. velocity it rms =< v 2 > 1//2 , the r.m.s. magnetic field b rms =< b 2 > 1 ^ 2 , the 
reduced cross-helicity p, the kinetic Reynolds number 1Z V = u rma L^jv (which is equal to the magnetic Reynolds number), the 
eddy turnover time r„i = L±/u rma and the Alfven time ta = Lu/Bq. All these quantities are also computed in the stationary 
phase for which turbulence is fully developed. Note that the value of ta for run la in the stationary phase is overestimated 
since b rm s > Bo.) We finally give the time t* from which the driving is stopped and tm the final time of the simulation. 
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C. 3D modes and 2D state 

In the presence of an external magnetic field Bo, it 
is convenient to describe the flow dynamics in terms of 
Alfven waves which propagate along Bo at frequency 
w(k) = k • Bo = fc||Bo- I n such a description, we may 
define the 3D modes, namely the spectral fluctuations at 
fcii =/= 0, for which the Alfven frequency is non-zero. The 
complementary part namely the spectral fluctuations at 
fen = defines the 2D state (or slow mode). Note that 
the 2D state can still be associated to waves since a local 
mean magnetic field bo may be defined in the fcn = plan 
along which local Afven waves propagate with frequency 
W2_o(k) = kj_ • b . 

To obtain the fields associated to the 3D modes and 
the 2D state a simple decomposition is performed from 
the Fourier space. For the 3D modes, the Elsasser, ve- 
locity and magnetic fields are defined as z ± (k±, ku > 0), 
u(fc_L,fcii > 0) and b(fej_,fe|i > 0) respectively. For the 
2D state, the velocity and magnetic fields are defined as 
u(fcj_,fc|| = 0) and b(A:j_,fc|i = 0). From these quantities, 
we define the Elsasser fields of the 2D state as 

z ± (k±,k\\ = 0) = u(kj_,k\\ = 0) ± b(kj_, k\\ = 0) . (14) 

In the rest of the paper, the quantities associated to the 
3D modes and 2D state will be noted with "w" and "2£>" 
indices respectively. Note that a different notation was 
used in |24| where shear- and pscudo- Alfven waves were 
defined from a toroidal/poloidal decomposition. 

III. TEMPORAL ANALYSIS 

In this Section, we study the temporal behavior of sev- 
eral global quantities to characterize the MHD flow dy- 
namics and the influence of the Bq intensity. In all the 
following Figures, time evolutions are shown for simula- 
tions la to Ilia and Illb (see Table H}. 

We first consider the evolutions of the total energy 

E(t)=E+(t) + E-(t), (15) 



where 

E ± {t) = \( [ 2±\k x ,k\{)dk x dk\\, (16) 

Z Jk n Jk ± 

the energy of the 3D modes 

E w (t) = E+(t) + E~(t) , (17) 

where 

Et{t) = \ [ [ 2 ±2 (fc±,fc||)dfc±&||, (18) 

1 Jk l{ >0 Jk ± 

and the energy of the 2D state 

E 2D (t) =E(t)-E w (t). (19) 

The temporal evolution of these energies are displayed 
in Fig. [T] For all simulations, we distinguish three parts 
in the evolution of the total energy E(t). The first part 
is characterized by an increase of the total energy which 
corresponds to a non stationary phase where the bal- 
ance between the external force and the dissipation is 
not reached. When it is reached a plateau is observed: 
it is the second phase. Finally, a third phase appears 
when the external force is suppressed at time t* (see Ta- 
ble HJ). Then, a decrease of the total energy is found 
during this decay phase. Note that the second phase is 
reached later when the mean magnetic field is stronger; 
it is compatible with a simple heuristic analysis in terms 
of time scales which shows that a transfer time in t^,/ta 
is indeed larger at larger Bq. 

The main difference observed between these three sim- 
ulations comes from the evolution of the energies E w (t) 
and E2D{t)- We remind that initially there is no energy 
in the 2D state, but as we can see, shortly after the be- 
ginning of the simulations the 2D state is excited. For 
the case Bq = 1, the energy increases mainly in the 3D 
modes whereas the energy of the 2D state reaches in the 
stationary phase approximatively 16% of the total en- 
ergy. A different behavior is found for the case Bq = 15 



4 



where a strong increase of the energy of the 2D state is 
measured which reaches in the stationary phase around 
2/3 of the total energy. In this case the energy of the 2D 
state is significantly larger than the one of the 3D modes 
which displays only a weak variation during the driving 
phase. The case Bo = 5 is an intermediate case where 
the energy of the 2D state stabilizes approximatively at 
1 /3 of the total energy. 
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FIG. 1. Temporal evolution of E(i) (solid), E 2D (t) (dash) 
and E m (t) (dot) for So = 1 (a), 5 (b) and 15 (c) (runs la to 
Ilia). Note that the vertical lines correspond to t* . 



It is important to note that during the first phase the 
increase of b rms leads to a decrease of the ratio Bo/b rms 
to around 0.4, 2.5 and 7.1 for runs la, Ila and Ilia 
respectively, with the characteristic timescale relations 
t~a ~ Tni for runs I — Ila and ta <§C t„j for runs Ilia (see 
Table |T|). This evolution has an impact on the nonlinear 
dynamics as we can see in run Ila: the energy of the 
3D modes has a weak variation until time t ~ 1.5 for 
which we have Bo/b rms ~ 3.3. After this period of time 
a stronger increase of the energy of the 3D modes is noted 
until a stationary phase is reached. The value Bo/b rrns ~ 
3.3 appears to be a threshold beyond which the energy 
in the 3D modes is roughly conserved. This analysis is 
confirmed by run Ilia where the energy of the 3D modes 
does not change significantly during the period of driving 
for which we always have Bo/b rms > 3.3. 



Figure [2] shows the time variation of the total energy, 
namely 



E(k h t) = E+(k h t) + E-(k h t) 



with 



fcii )dfej 



(20) 



(21) 



where fen is fixed to 0, 1, 2 and 4. We first note that the 
energy in plan feii = 1, 2 does not change very much at 
Bo = 15 whereas it does change for other values of Bo- 
This behavior is surprising since the external force acts in 
particular on plans fen = 1 and 2. This situation contrasts 
with the energy in the fey =0 plan where the total en- 
ergy exhibits a strong increase for Bo = 15. For the plan 
fey = 4 we note that the energy reaches only a small value 
when Bo = 15. This observation demonstrates that the 
2D state pumps efficiently the energy injected into the 
system at small parallel wavenumbcrs when a strong ex- 
ternal magnetic field is applied. The second interesting 
observation in Fig. [5] is the oscillations found in the plan 
feii = 4 (and also for higher values of k\\ ; not shown) at 
the very beginning of runs Ila and Ilia. During this pe- 
riod of time which extends approximately up to t = 1.5 
for Bo = 15 and t — 1 for Bo — 5, the ratio Bo/b rms 
is larger than 10 and 4 for respectively Bo = 15 and 5. 
These oscillations are apparently correlated to the inten- 
sity of the mean magnetic field with a higher frequency 
at stronger Bo/b rms . The last remark is about the dy- 
namics at time t > t*. For all simulations and all energies 
we note generally a decrease in time - sometimes sharp 
like in fen =1,2 plans for Bo = 1 - which corresponds 
to a decay phase. Note that for Bq = 15 the decay is 
characterized initially by small oscillations which are in 
phase opposition between the energies in plan fey = and 
1. This feature illustrates the energy exchange between 
the two first fey plans. 

Figure [3] presents the Alfven ratios between the kinetic 
and magnetic energies for the 3D modes 



the 2D state 



and the total flow 



rib® 



r A {t) 



E v 2D (t) 
E b 2D it) ' 

E v (t) 
E~Hi) ' 



(22) 



(23) 



(24) 



where the definitions (TTB]), ([T5]) and ([TO]) have been di- 
rectly applied to the velocity and magnetic fields. The 
Alfven ratio allows us to measure the prevalence of Alfven 
wave fluctuations: for example, in the wave turbulence 
regime one can demonstrate at the level of the kinemat- 
ics an equipartition between the kinetic and magnetic 
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FIG. 2. Temporal evolution of E(k\\,t) for B = 1 (dot), 5 (dash) and 15 (solid) (runs la to Ilia) at fcy = (a), 1 (b), 2 (c) 
and 4 (d). Note the use of semi-logarithmic coordinates for fey = 4 (bottom right). 



energies [l4j . Therefore, a departure from unity suggests 
the presence of non Alfvenic fluctuations [24[ . In simu- 
lations Ila and Ilia, the 3D modes are in equipartition 
with ~ 1. However, we may note a difference be- 
tween t < t* and t > t*: a lack of oscillations is found 
during the driving phase whereas significant oscillations 
happen during the decay phase which are stronger for 
Bq = 15. The 2D state evolves quiet differently with an 
Alfvcn ratio r^ D significantly smaller than the unity for 
runs la to Ilia which means that the 2D state is magnet- 
ically dominated. We also note that r^ D seems not to be 
strongly affected by the value of the mean field Bq with 
first an increase (from the initial value at 0.4) and then a 
slight decrease. The behavior of the Alfven ratio for the 
2D state is therefore quiet similar to what was found in 
the pure decay regime (24|. The time variation of r A is 
smooth with initially only a slight decrease after which 
the value stabilizes around 0.8 for Bo = 5 and 15. When 
the driving is suppressed a slight decrease with small os- 
cillations is found as expected since it is the addition of 
the effects of r A and r^u- For run la where Bq — 1 a 
different behavior is observed since all Alfvcn ratios are 
less than unity and no Afvenic fluctuations are detected. 
This behavior is partly explained by the increase of b rms 
which becomes larger than Bq and prevent any dynami- 
cal effect of the mean magnetic field on the flow. 

Figure S] shows the reduced cross- hclicity p(t) between 
the velocity and magnetic fields which was already de- 
fined by relation (jlOp . In addition, we may define the 



reduced cross-helicity for the 3D modes and the 2D state 
as respectively 



and 



where 



Pw(t) 



P2D(t) = 



E+(t)-E~(t) 
E+(t)+E^{t) ' 



E+ D {t) + E- D (t) 



Efoit) = \ I z ±2 (fc±,fc|| = 0)dkj 
1 Jk 



(25) 



(26) 



(27) 



We remind that the cross-helicity is a measure of the 
relative amount of Alfven wave packets which propagate 
in opposite directions along the uniform magnetic field 
Bo- It is generally found that MHD turbulence evolves 
towards a state of maximal cross-helicity with an align- 
ment or an anti-alignment of v and b according to the 
sign of the cross-helicity 0- This definition applies to 
(f2"5"|) which is a refined definition of the cross-helicity 
where the non propagating part is discarded. With the 
definition (|2"fj)) the situation is different. If we assume 
isotropy in the feii = plan then p2D can be seen as a 
local measure of the cross-helicity since we can always 
define a local mean magnetic field bo along which 2D 
Alfven fluctuations propagate. In Fig. @] we see that the 
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FIG. 3. Temporal evolution of Alfven ratios r A (t) (solid), 
r-2 D (t) (dash) and r„(t) (dot) for Bo = 1 (a), 5 (b) and 15 (c) 
(runs la to Ilia). 
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FIG. 4. Temporal evolution of p(t) (solid), P2D(t) (dash) and 
p w (t) (dot) for Bo = 1 (a), 5 (b) and 15 (c) (runs la to Ilia). 



reduced cross-helicity of the 3D modes is initially null 
and slightly evolves during the driving phase for Bq = 1 ; 
for B = 5 only small oscillations are detected whereas 
for Bo = 15 it is almost constant. One has to wait the 
decay phase to see a significant modification of p w . 

The situation is different for the reduced cross-helicity 
of the 2D state which experiences a variation whatever 
the intensity of the mean field Bq is. Note that this vari- 
ation is moderate since it never exceeds in absolute value 
0.35. Note also that our initial condition leads to a neg- 
ative value for p2D since we impose a condition only on 
the global cross-helicity, i.e. p(t = 0) = [p2D cannot 
be defined initially since the 2D state has no energy). 
In other words, we generate initially more z^ D than z^ D 
fluctuations. For Bq = 5 and 15 the time at which the 
stationary phase is reached are t ~ 4 and 5 respectively; 
at these times, pio indicates around 30% more z^ D than 
Z2 D - Clearly, we see that the evolution of p(t) is mainly 
affected by the evolution of p2D(t)- Therefore, it seems 
to be important in this problem to make the distinction 
between the cross-helicity of the 3D modes and the 2D 
state to evaluate if whether or not we have the domina- 
tion of one type of Alfven wave packet propagating along 
the external magnetic field Bq. Another interesting be- 



havior is the variation of p2D during the decay phase: we 
observe a weaker variation for Bq = 5 and 15 with in the 
latter case an almost constant reduced cross-helicity. 

IV. SPECTRAL ANALYSIS 

In this Section, we study the spectral behavior of 
driven MHD turbulence. We start from the bidimen- 
sional Elsasser energy spectra, E ± (k±, ku), and define re- 
spectively the corresponding unidimensional spectra for 
the 3D modes 

Et{k±)= I E ± (k ± ,k ll )dk ll , (28) 
A-n>o 

and the 2D state 

E± D (k ± ) = E ± (k±,k ll =0). (29) 

Then, the addition of both defines the usual unidimen- 
sional Elsasser energy spectrum 

E ± (k ± )=[ E ± (k ± ,k ll )dk {l . (30) 

The comparison will be made between all runs, i.e. runs 
la to Ilia and run Illb. 
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FIG. 5. Energy spectra E + (k±) (solid), E^ D (k±) (dash) and E+(k±) (dot) for B = 1 (top), S = 5 (middle) and B = 15 
(bottom) which corresponds to runs la, Ila and Ilia respectively. Spectra are displayed at times t < t* (left) and t > t* (right) 
which corresponds to t = 5 and 7.2 for Bo = 1, t = 5.1 and 9.2 for Bo — 5, t = 5.8 and 6.2 for _Bo = 15. 



A. Runs la to Ilia 

Figure [5] displays energy spectra for the Elsasser 
field E+(k±), the 3D modes E+(k±) and the 2D state 
E^ D {k±), at different Bo intensity and with the same 
initial condition (runs la to Ilia). The comparison is 
made between two times, t < t* and t > t* , where t* 
is the time from which the external force is suppressed. 
Shortly after the beginning of the simulations we observe 
that the energy spectrum i?2z>(fcx) loads up. For Bq = 1 
this spectrum remains for all times about one order of 
magnitude smaller than E+(k±). The situation is dif- 
ferent for stronger magnetic intensity like Bq = 15. In 
this case, we see that £ , ^ D (/c^) may dominate E+(k±) at 
small perpendicular wavenumbers, i.e. k± < 10. Then, 
the spectrum E + (ki_) is formed at large scales mainly 
by the 2D state and at small scales by the 3D modes. 
Run Ila is an intermediate case where the 2D state be- 
comes dominant in a very narrow range of perpendicular 
wavenumbers. The second interesting comment is about 
the break observed in some 3D modes spectrum just af- 



ter the domain of excitation which is limited in particular 
to k± e [1,2] (see relation (J7J). This break appears for 
strong Bo and during the driving phase. Indeed, shortly 
after time t* an extended power law energy spectrum is 
formed. At this resolution the driving prevents the pos- 
sibility to measure any power law for the energy spec- 
trum of the 3D modes; it is believed that it could also 
artificially modify the power law of the Elsasser energy 
spectra. Note that a similar break was already observed 
in driven MHD turbulence [36] - Finally, at time larger 
than t* the system evolves freely and the slope at the low- 
est k± is modified slightly for strong Bq and strongly for 
Bq = 1 with a flattening for the energy spectrum of the 
3D modes: the slope goes from k\ initially to roughly k\ 
after t*. Note that E~(k±), E~(k±) and E^ D (k±) (not 
shown) exhibit a similar behavior. 

The bidimensional Elsasser energy spectra E + (k±, kt\) 
are plotted in Fig. [6] at fixed k\\ . We choose the values 
fcii = to 5 and also fc|| = 15. As in Fig. [5] a comparison 
is made between two times, t < t* and t > t*. The most 
remarkable feature is that the spectra at fixed fey are 




FIG. 6. Energy spectra E + (k±, k\\) at fixed k\\ and for Bq — 1 (top), £?o = 5 (middle) and Bo = 15 (bottom) which correspond 
to runs la, Ila and Ilia respectively. Spectra are displayed at times t < t* (left) and t > t* (right) which are the same as in 
Fig- El The values chosen for fey are 1, 2, 3, 4, 5, 15. For comparison we also plot the spectra of the 2D state (green color online, 
grey line) for which k« = and of the 3D modes. 



characterized by a pinning effect with a convergence of 
the spectra at large k±. The case Bq = 15 is an exception 
where the pinning is observed only for close values of fen 
as we can see with a significant difference for fcii = 15. 
Note that the spectrum of the 2D state has a slightly 
different behavior since it does not converge towards the 
other spectra at large k±. The second observation is 
that at times larger than t* power laws are detected for 
spectra with fen = to 5, with a reduction of the inertial 
range while the power law index seems unchanged. The 
inertial ranges are then shifted to larger k± . 

Figure [7] presents unidimensional Elsasser energy spec- 
tra E ± (k±) defined by relation (|30p. These spectra are 
compensated by different power laws k™ with m = 3/2, 
5/3 and 2 which correspond to the predictions for strong 
and weak turbulence. Only runs Ila and Ilia are shown 
for which respectively Bo/b rrns = 2.5 and 7.1 (run la 
behaves similarly to run Ila). We see that run Ila fits 

well with power laws like k ± 3 ^ 2 or k ± 5 ^ 3 whereas run 
Ilia is close to fcj 2 . These results are in relatively good 



agreement with the characteristic timescales given in Ta- 
ble U which satisfy relations ta ~ t„; for run Ila and 
ta <C r n i for run Ilia. We also see that for run Ilia the 
inertial range is significantly larger for E + than E~ . We 
will come back to this point below and show that the 2D 
state may explain such a feature. 



B. Runs Ilia and 111b 

In this Section a comparison is made between runs 
Ilia and Illb for which Bq = 15. The difference be- 
tween these runs comes from the initial conditions and 
the external force (see Table U and Sec. Ill Bli . 

Figure E] shows the compensated energy spectra 
E ± (k±) for run Illb. As in Fig. [7] spectra are com- 
pensated by different power laws fc™ with m = 3/2, 5/3 
and 2. We observe the same characteristics as for run 
Ilia with in particular timescales compatible with rela- 
tion ta <C r n i (sec Table QJ. The only one difference is 
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FIG. 7. Energy spectra E ± (k±) compensated by k\ (solid), k 3 / 2 (dash) and k 5 / 3 (dot) at times t < t* (left), t > t* (right), for 
runs Ila and Ilia. Note that the times are the same as in Fig. [5] 
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FIG. 8. Energy spectra E (k±) (run Illb) for times t = t* =5.1 (left) and t = 5.2 (right) compensated by k ± (solid), 
(dash) and k 5 / 3 (dot). 
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the size of the inertial ranges which are about the same 
for both spectra. 

In order to explain the previous observations and the 
difference between cases Ilia and Illb we shall analyze 
the ratios E^/E^d- We first remind that in the regime 
of wave turbulence the wavevectors (k, p, q) satisfy the 
resonant triadic interactions k = p + q and for example 
w(k) = w(p) — w(q) which leads to the resonance con- 
dition g|| = and to the prediction that only a spectral 
transfer transverse to Bo happens. It also means that 
we always have an interaction between an Alfven wave 
packet of one type of polarity and the fluctuations of 
the 2D state with the opposite polarity. Numerical sim- 
ulations have confirmed that the resonant interactions 
become dominant for large B$/brms [Hi- Therefore, the 
amount of energy in the 2D state is an important pa- 



rameter for the nonlinear dynamics. For example, the 
absence of energy inside the 2D state cannot lead to the 
development of a wave turbulence regime. In Fig. [9] the 
ratios E^ / E^jj between the energies of the 3D modes and 
the 2D state with opposite polarities are given. For run 
Illb, we see that these ratios are almost the same and 
reach roughly a stationary state for times t £ [4, 6] with 
values close to 1.5. In this regime we have E+ D ~ 0.59E- 
and E^ D ~ 0.66.E+ which lead to similar spectra E ± (k±) 
(see Fig. |8]) with inertial ranges of the same size. For run 
Ilia the ratios are significantly different: in the station- 
ary phase for which t £ [4, 6] we have E^ D ~ 0.33i?~ and 
~ ^w- in practice, the energy in the 2D state with 
a positive polarity is reduced which leads to a weakening 
of the nonlinear dynamics and a reduction of the inertial 
range for E~(k±) as we can see in Fig. [71 Note that sim- 
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FIG. 9. Temporal evolution of the ratios E^/E 2D for runs 
Ilia (top) and Illb (bottom). 



est k±, the initial power law in k\ applied to spectra 
fcii = 1, 2 is slightly modified and the scaling for E^(k±) 
is unchanged. Note also that the break previously re- 
ported for times t < t* is also observed (not shown). 
Finally, as we have already seen in Fig [5] (run Ilia) , for 
times t > t* an inertial range appears for spectra E^(k±) 
and E (k±,k\\) with k\\ — 1,2,3,4 which seem roughly 

better fitted by k ± 3 ^ 2 or k ± 5 ^ 3 than kj 2 . 
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ilarly to run Ilia (see Fig[I) the energy contained in the 
2D state for run Illb represents 2/3 of the total energy. 




1 10 100 

ki 



FIG. 10. Energy spectra E + (k±, k\\) (top) and E~(k±, fe||) 
(bottom) at fixed ku (run Illb) and at time t — 5.5 > t*. 
The values chosen for ku are 1,2,3,4 (same notation as in 
Fig. [6]). For comparison we also plot the spectra of the 2D 
state (solid) and of the 3D modes (bold solid). 

Figure [10] displays spectra E^(kx), E 2D {k±) and 
E (fcj_,fc||) for /c|| = 1,2,3,4. Globally we observe the 
same features as for run Ilia with a stronger domination 
of the 2D state at large scales. Note that at the small- 



FIG. 11. Temporal evolution of Alfven ratios r A (t) (solid), 
r^oit) (dash) and r£(t) (dot) for run Illb. 

Figure [TT] shows eventually the time variations of the 
Alfven ratios r A (t), r^uit) and r^(t). The main differ- 
ence between run Ilia (see Fig. [3]) and run Illb resides 
first in the initial values which are larger for the latter 
run. The evolution is also different with a variation of 
all quantities during the driving phase. It is only during 
the decay phase that we recover the same behavior as in 
run Ilia with alfvenic fluctuations and an equipartition 
for the 3D modes. In this decay phase, a large devia- 
tion of r^(t) is observed immediately after t* in order 
to reach the equipartition which clearly means that it 
is the natural state of the system (in Ilia the situation 
was not totally clear since before t* we were already at 
equipartition) . 

V. SUMMARY AND CONCLUSION 

In this paper, we present a set of 3D direct numeri- 
cal simulations of incompressible driven MHD turbulence 
under the influence of a uniform magnetic field Bo- In 
particular, the temporal and spectral properties of the 
2D state (or slow mode) are investigated. We show that 
if initially the energy contained in the 2D state is set 
to zero it becomes shortly non negligible in particular 
when the intensity of Bq is strong. For our larger Bq 
intensity (Bq = 15) this energy saturates in the station- 
ary phase around 2/3 of the total energy whereas the 
energy of the 3D modes remains roughly constant. We 
also observed that the ratio Bq /b rms saturates around 7 
for both simulations with Bq — 15 with initially a ratio 
equal to 15 (run Ilia) and to 19 (run Illb). In all situa- 
tions, the magnetic energy dominates the kinetic energy 
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but it is shown that at large Bo/b rms and in the decay 
phase the natural state for the 3D modes is the equipar- 
tition whereas the 2D state is magnetically dominated. 
It is interesting to note that a theoretical model in terms 
of condensate has been recently proposed to explain the 
spontaneous generation of a residual energy E v — E b at 
small fen [I?}. This model is based on the breakdown of 
the mirror-symmetry in unbalanced wave turbulence. 

From a spectral point of view, when the B intensity 
is strong enough the k j_ -energy spectra are mainly com- 
posed at large scales by the 2D state and at small scales 
by the 3D modes. This situation is similar to rotating 
turbulence for neutral fluids [13, HH where in particu- 
lar the nonlinear transfers are reduced along the rotating 
rate (39|. However, a detailed analysis of the temporal 
evolution of the 2D state energy spectra shows that a 
direct cascade happens whereas for rotating turbulence 
an inverse cascade is generally evoked [4^] . The same re- 
mark holds for energy spectra at fixed fcii > where we 
observe additionally a pinning effect at large k± . Accord- 
ing to the value of So/6 rms scalings close to k ± 3 ^ 2 - k ± °^ 3 
or kj 2 are found which are in agreement with different 
predictions for strong and wave turbulence 0, [^, [2 . 

The external force seems to be an important param- 



eter for the dynamics. For example in |25[ a change of 
spectral slope was reported for the k± -energy spectrum 
when the intensity of -Bo is modified. In this work the 
external force was applied at ku > and kj_ = 1,2. In 
[lH a forcing which kept the ratio of fluctuations to mean 
field approximately constant was implemented by freez- 
ing modes k < 2 and a spectral slope close k 3 ^ 2 was 
reported for Bq = 5. The fact that the 2D state was 
imposed by the forcing has certainly an impact on the 
nonlinear dynamics since it prevents the natural growth 
of the 2D state energy at small k±. We remind that 
the 2D state is essential at large Bo/b rrns since it mainly 
drives the nonlinear dynamics. In particular, we observe 
that a reduction of the 2D state energy of one type of 
polarity leads to a decrease of the inertial range of the 
3D modes energy spectrum with the opposite polarity. 
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